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Abstract 

Charmonium suppression in hot and dense nuclear matter has been argued to be a signature 
for the production of the quark gluon plasma (QGP). In order to search for this effect in heavy 
ion collisions one must have a clear understanding of all the factors that can contribute to such a 
suppression. These may include shadowing of the partons in a nuclear environment, breakup of 
a correlated c-c pair as it traverses the nuclear fragment, suppression of feed-down from higher 
mass states as well as other initial state interactions. In order to disentangle these effects one 
must measure charmonium production rates in both proton+proton (p+p) and proton+nucleus 
(p+A) collisions. The p+p collisions serve as a baseline for searching for suppression compared 
to binary scaling predictions, allow one to quantify the amount of feed-down from higher states 
as well as serve as a tool to distinguish between different theoretical calculations for charmonium 
production mechanisms. In order to quantify nuclear effects it is also necessary to study char- 
monium production in p+A collisions where the temperature and density of the system are low 
compared to a heavy ion collision. These measurements allow one to determine the influence of 
nuclear shadowing and breakup in "cold" nuclear matter which can be extrapolated to heavy ion 
collisions in order to determine the amount anomalous suppression. Of course, extrapolations 
that rely on a model based technique depend heavily on the assumption of a production mecha- 
nism, a fact that reinforces the importance of the p+p measurements. The PHENIX and STAR 
experiments at Brookhaven National Laboratory have measured charmonium production in p+p, 
d+Au, Au+Au and Cu+Cu collisions at y/s^N - 200 GeV for both forward and mid rapidities. 
I will present a review of the latest measurements from both experiments with an emphasis on 
what we have and can still learned from them about charmonium production and suppression 
with these experimental apparatuses. 



1. WhytheJ/t^ 

The heavy nature of charmonium (cc) allows one to apply potential models in non-relativistic 
quantum mechanics to calculate the mesons binding radius. Originally it was predicted that the 
modification of the heavy quark pairs potential in the hot dense matter created in heavy ion 
collisions would cause the pair to become uncorrelated due to color charge screening. This 
modification of the pair potential via a Debye mass term leads to charmonium suppression 1 1 1 
when compared to a binary collision scaled p+p reference. Due to the different binding energies 
for the different charmonium states one could gain access to the temperature of the medium. At 
RHIC the suppression of the lowest energy charmonium state, the J/iff meson, has been measured 
in Au+Au and Cu+Cu collisions in V^aw =200GeV collisions. 
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These measurements are then compared to the invariant yield measured in a baseline p+p 
collision at the same center of mass energy. It is assumed that any modification due to the 
medium will show up as a deviation from the prediction of the binary collision scaled reference 
data. This is of course true when no modifications due to normal nuclear matter are present. 
However in the case of the J /iff we know this is not the case from lower energy measurements 
made in p+A collisions at SPS [2] and FNAL 0. 

Therefore, to interpret these data one must remove any effects that occur in normal density 
cold nuclear matter (CNM). One such effect is the modification of parton distribution functions 
in a nuclear environment |j4j O . Another is Cronin enhancement that leads to a hardening of the 
transverse momentum spectrum of collision products due to multiple scattering. To this end the 
PHENIX experiment has also measured the nuclear modification present in d+Au collisions in 
a/* =200GeV collision. The d+Au data is used to extrapolate within a Glauber [6| based data 
driven model to the Au+Au collision case to predict the suppression that would result from CNM 
effects and search for an anomalous suppression of J /if/ mesons in the sQGP. 

2. J /iff Production 

The production mechanism for charmonium has not been well understood theoretically for 
nearly 20 years. The magnitude of the pj spectrum measured at the Tevatron Q|8| was under- 
predicted by an order of magnitude by the color singlet (CS) model [5J. This led to the proposal 
of the color octet (CO) model flOl fTTI wherein the pre-charmonium charm quark pair become 
correlated in a color charged state and must color neutralize via soft gluon emission. While 
the CO model had some success in describing the kinematic spectra it also predicted a large 
transverse polarization at intermediate to large pj IfTUl [TP . Another undesirable aspect of the 
CO is its reliance on overlap matrices between the color charged pre-hadron and final state J/if/, 
which are nearly free parameters in the model fl2l . 

Measuring the p+p spectra and polarization at RHIC is a tool for distinguishing between 
charmonium production mechanisms. PHENIX has measured the invariant yield of J/if/s, in p+p 
collisions at V* — 200 GeV, over a wide range in transverse momentum at forward (1.2 < \y\ < 
2.2) and mid rapidity (|y| < 0.35) (Figure [T] top) as well as the J/iff polarization at mid rapidity 
(Figure [T] bottom). The quality of these data make them natural metrics for testing new models 
of charmonium production such as the four-point modified CS model proposed by Lansberg and 
Haberzettl lfl3l . The mid-rapidity data agrees well with the polarization predictions. At forward 
rapidity there is a two sigma difference between the measurement and prediction. 

Another, quite recently measured, observable that may help establish the production mech- 
anism for the J /if/ is azimuthal correlations with hadrons. One could imagine that the spatial 
correlation between the }/ifr and the remainder of the hadrons forming the jet that results from 
the hard initial collision may be very different depending on the production mechanism. STAR 
has recently released [16| azimuthal correlations between }/if/s and hadrons(Figure|2]l measured 
in y/siw = 200 GeV p+p collisions. The J/if/ trigger particle is required to have pj > 5 GeV/c 
and the associated hadrons to have pj > 0.5 GeV/c (see Section 13). This new observable may 
prove very fruitful in distinguishing between different production models in the future. One can 
imagine handing a MC generator (i.e. PYTHIA) the QCD calculation for different models and 
extracting the resulting correlation of prompt ]/if/ mesons and hadrons within experimental cuts 
and then comparing these to data. 

Understanding the production mechanism establishes the map between the measured kine- 
matics of the 3/if/ (y, pj) and the kinematics of the partons (xu X2) at the interaction vertex. Dif- 
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Figure 1 : (color online) Invariant yield of ]/if/ measured with the PHENIX spectrometer versus pj (top left) and y (top 
right). The theoretical curves are the four point modified CS model [141 and NRCD 1 15 ]. J/t/< polarization versus pj 
(bottom left) measured at mid and forward rapidity with the PHENIX spectrometer. The dashed lines are the predictions 
from the four point modified CS model 1 14 |. 



ferent production models provide these mappings and can result in different conclusions about 
the magnitude of the cold nuclear matter effects ifTTl discussed in Section [4] The J /iff polariza- 
tion and azimuthal correlation measurements provide us with another lever arm to distinguish 
between these models. 



3. Feed Down 

One of the contributors to the J/t/> spectrum is expected to come from the feed down of higher 
charmonium states as well as B-meson decays. The tf/\ifr(2S)), \c and B mesons all have a decay 
mode to J/i/z+X. Considering the predictions from the lattice for the dis-association temperatures 
of these states at < T c Ifl8ll (where T c » 170 MeV |19|) and the most recent interpretations of 
the temperatures reached in heavy ion collisions at RHIC ~ 1.57c [20] one would expect the 
depletion of at least the;tv and the t/>' which would in turn lead to a J/t/> suppression beyond that 
in cold nuclear matter collisions. Of course this depletion must also be accounted for when one 
attempts full accounting of ]/ifr suppression in heavy ion collisions. 

To this end, the PHENIX collaboration has measured the feed down fraction (R) of the i//' 
and set a 90% confidence level upper limit for the \c to J/i/> in p+p collisions at 200 GeV via 



3 



the di -electron channel at central rapidity 11211 . The results are in agreement with theoretical 
predictions Il22l as well as the world average [23 1 (see table in [24]). 

As mentioned above the STAR experi- 
ment has also measured the azimuthal cor- 
relation of J /iff and hadrons (Fig. [2}. Using 
this correlation one can extract the feed-down 
fraction of J /iff arising from B meson decays 
if one knows the azimuthal shape of the B- 
decay and prompt J /iff correlations. In order 
to extract this quantity STAR parameterizes 
the shape as C(A<p) = x * C p (A(f>) + (1 - x) * 
Cb(A<P), where C p is the correlation function 
from prompt J /\fj mesons and Cb comes from 
the B decay. The prompt correlation function 
is then taken from a PYTHIA simulation that 
has been tuned (i.e.both color octet and sin- 
glet mechanisms mixed) to match real data 
transverse momentum distributions. While this may be a rather drastic assumption, it is nec- 
essary to extract the B feed down fraction in this manner. One would naturally argue that an 
assumption of this nature would inherently introduce a large systematic uncertainty. This is es- 
pecially true if one considers that the production mechanism (as discussed above) is still not 
known and that PYTHIA simulations do not reproduce either the cross section or polarizations 
measured in data. However, this extraction yields a feed-down fraction of 13 ± 5 % lfl6ll . 

4. Cold Nuclear Matter Effects 

In order to interpret suppression in heavy ion data one must remove any effects that occur in 
normal density cold nuclear matter (CNM). One such effect is the modification of parton distri- 
bution functions in a nuclear environment |4]|5). Another is Cronin enhancement that leads to a 
hardening of the transverse momentum spectrum of collision products due to multiple scattering. 
The PHENIX experiment has also measured the nuclear modification present in deuteron+gold 
(d+Au) collisions in -\fs — 200 GeV collision. 

During the 2008 RHIC deuteron-gold (d+Au) run PHENIX recorded a factor of 30 greater 
integrated luminosity than the previous run in 2003. The PHENIX collaboration has not yet 
released nuclear suppression factors for d+Au collisions compared to the Run-3 p+p baseline. 
However, the data were analyzed to measure the central to peripheral ratio R CP in these collisions 
(Eq.[T]i. In order to calculate the central to peripheral ratio Rqp the invariant yield of J /iff mesons 
from d+Au central and peripheral collisions must be measured. Figure [3] shows the nuclear 
modification as a function of rapidity (y) in a given centrality bin (/) which defines the average 
number of binary collisions (N"'"). 

These data have sparked much interest in the community (for details see the 2009 ECT and 
INT quarkonia workshops) as it is clear that a rapidity independent breakup cross section com- 
bined with nuclear shadowing cannot match the shape of the data. This has lead to the conclusion 
that there may be some physics missing in the models. PHENIX is currently working to under- 
stand a normalization effect, between different RHIC runs, for the nuclear suppression factor 
(RdAu) an d these results are expected to be made available later this year. 




4 



PHENIX preliminary 




-RHIC 2008 

d+Au\Js^ = 200 GeV 
"±11% Global Scale Uncertainty 

40-60% 




±5.05% Global Scale Uncertainty 

0-20% 



ranp 



|y|<0.35 syst =±12% 

" ' ' global 

|y|e[1.2,2.2] syst =+7% 



ffl 



i 



50 100 150 200 250 300 350 400 



• STAR Cu*Cu 0-60% 

• PHENIX Cu+Cu 0-94% 



Figure 3: (color online) Central to peripheral ratio (Rc p ) measured by the PHENIX collaboration using the RHIC 2008 
d+Au data set (left). The theoretical curves represent the prediction for Rcp using the EKS 1 25 1 nPDF parameterization 
and different constant-with-rapidity breakup cross sections. Published PHENIX Au+Au nuclear modification factor |26| 
showing a larger suppression at forward rapidity than central (top right). Comparison of the PHENIX and STAR results 
for nuclear modification in Cu+Cu at V* = 200 GeV (bottom right). 
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5. What's Hot in Hot nuclear matter 

In Figure [3] (top right) we show the nuclear suppression factor for Au+Au collisions at 
= 200 GeV, as a function of centrality, measured at forward and mid rapidity by the PHENIX 
collaboration at RHIC (26). One striking feature of this comparison is the similarity between 
the suppression patterns at mid rapidity between the RHIC (|y| < 0.35) and SPS (0 < y < 1) 
data [27), despite the difference in center of mass energies of the two measurements (yfs — 200 
GeV and yfs ~ 20 GeV respectively). This engenders the question "Why are the results so 
similar?" for measurements probing different rapidity regions, shadowing regimes and energy 
densities. It is clear from the PHENIX data that the suppression at forward rapidity is greater 
than at mid rapidity. This is a challenge to local density based suppression models and is one 
piece of evidence that supports the idea of regeneration discussed below, where the J /iff yield is 
enhanced due to close proximity in phase space of uncorrelated pairs. 
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Figure 4: (color online) Elliptic flow (V2) versus pj measured with the PHENIX detector at mid and forward rapidity 
compared to various model predictions 1 33 34 35 36 371 ■ 



We also show (Figure [3] (bottom right)) the high pj measurements of Raa made by the 
PHENIX and STAR experiments in Cu-Cu at ^Js^ = 200 GeV. One should note that the 
PHENIX measurement is for minimum bias collisions while the STAR measurement is for cen- 
trality 0-60%. Given this small difference one can still compare the measurements. It is clear 
that the STAR data favors a reduction of the suppression at high pj (with large systematic and 
statistical uncertainty), while the PHENIX data favors a nearly constant suppression from mid to 
high p T . 



6. A glimpse of the future. 

Regeneration (also called coalescence) |28l |29l is a process whereby un-correlated cc pairs 
coalesce to form charmonium states, an effect that could be enhanced due to the close phase 
space proximity of partons in a heavy ion collision. This process would increase the yield of 3/if/ 
in heavy ion collisions and must be accounted for when interpreting the overall magnitude of 
suppression. 

It has been suggested that one metric for determining the magnitude and/or existence of 
this effect would be to measure the elliptic flow (v'2) of the J /if/ in heavy ion collisions. The 
elliptic flow of the recombined J/if/ would be related to the flow of uncorrected charm 1 30] in the 
medium and could be very different from that of iif/ formed early on in the fireball's evolution. 
This experimental signature can be compared to the V2 from open charm and models that contain 
the regeneration mechanism. The PHENIX experiment has made the first measurement of J /if/ 
elliptic flow at forward and mid rapidity QTI for RHIC energies. The result does not have the 
statistical precision to distinguish between models and the comparison to the open charm elliptic 
flow l32l is not very enlightening (Figure [4j. A long, high-luminosity Au+Au run at RHIC will 
be necessary for this metric to be improved. 

In addition STAR has measured a clear T peak in the d+Au data taken in 2008 Figure [5] 
(right). This coupled with the previous p+p measurement [38] can be used to calculate RdAu f° r 
the T. The value is 0.98 + 0.32 (stat.) ± 0.28 (syst,). This measurement suffers from very large 
uncertainties which engenders the question of whether or not the RHIC program for higher mass 
charmonium states would benefit from a long d+Au run in the future. The current measurement 
would allow for suppression in cold nuclear matter of up to 30% within one standard deviation of 
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the statistical error alone. A long d+Au run will be required if one wants to precisely characterize 
any cold nuclear matter effects that would be present for the T state. 
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Figure 5: (color online) Nuclear modification probability distribution from PHENIX for T suppression in Au+Au 
collisions (left). Di-electron invariant mass spectrum from STAR using the Run-8 d+Au data that shows a prominent T 
peak (right). 



PHENIX has also released a new measurement of T suppression in Au+Au collisions. Due 
to a lack of statistics the value is presented as an upper limit on the amount of suppression seen 
in a HI collision. An upper limit of 0.64, at a 90% confidence level, is calculated by coupling the 
Poisson probability distributions for the numerator (Au+Au) and denominator (p+p) to arrive at 
a probability distribution for RauAu Figure. [5] (left). Interpreting this result is difficult due to the 
lack of information regarding cold nuclear matter effects for the Y. However, note that if the 
T(2s) and T(3s) states are melted in the medium a suppression of 30% would be expected in HI 
collisions. 
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